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Very thick films of BiCr03 have been grown on a SrTi03 substrate, maintaining a tetragonal lattice up to 
thicknesses of 14,000A. Assuming we can treat films of this thickness as bulk crystals, we first calculated 
the experimentally undetermined atomic positions within the unit cell with the measured lattice constant 
of the film, then relaxed the lattice constants along with the atomic positions. All the calculations result 
in an antiferroelectric crystal with the Pnma space group with noncollinear Bi displacements. We find 
G-type antiferromagnetism with a spin-orbit induced weak ferromagnetic component, however, the weak 
ferromagnetic component can cancel if the antiferromagnetic spins are oriented along a particular one of the 
three 2-fold rotation axes. 

PACS numbers: 77.55.Nv, 75.50.Ee, 75.70.Tj, 75.85.+t 



I. INTRODUCTION 

Ferroelectric antiferromagnetic crystals often have 
a weak spin-orbit induced ferromagnetic component. 
Fennie^ has recently shown that for a particular class of 
these materials, a reversal of the ferroelectric polarization 
results in a reversal of the ferromagnetic polarization. 
This gives these materials potential device applications 
as well as scientific interest. Unfortunately, the ferromag- 
netic moment is probably too small to be of practical use 
and the Currie temperatures are usually well below room 
temperature^. Dzyaloshinsky^ postulated that weak fer- 
romagnetism could occur only if the ferromagnetic state 
maintained the magnetic (color) group symmetry of the 
antiferromagnetic state. For a crystal with a trigonal 
axis, he found that D • Si x S2 was an allowed term in 
the Hamiltonian which favored spin canting, where Si 
and S2 are neighboring spins and D is a vector along 
the single trigonal axis. Moriya^ then found a formula 
for D. We will show how the Dzyaloshinsky-Moria (D-M) 
theory manifests itself in a crystal with three 2-fold sym- 
metry axes, i.e. how is the 2- fold-axis along which the 
D vector points to determined? Does that axis depend 
upon the orientation of the spins prior to the inclusion of 
the spin-orbit coupling(SOC)? 

BiCr03(BCO) was first synthesized over forty years 
ago by Sugawara et alP* Their finding of a triclinic low 
temperature crystal structure is, however, in conflict with 
some, but not all, more recent monoclinic determina- 
tions. There is a structural phase transition at about 
430 K, above which it agreecP*^ that BCO crystallizes in 
an orthorhombic structure with the Pnma space group 
(no. 62). Below 430 K x-ray powder diffraction patterns 
(XPD) were interpreted^ to yield a C2(space group no. 5) 
monoclinic structure. The lack of inversion symmetry 
implies a ferroelectric crystal. While, on the other hand, 
two different neutron powder diffraction (NPD)3S1 exper- 
iments led to the conclusion that, although the lattice pa- 
rameters were almost identical to XPD datcP, the space 
group is the inversion symmetric C2/c(no.l5). This cou- 



pled with the observed positions of the Bi atoms lead 
to the conclusion that BCO is antiferroelectric. G type 
antiferromagnetism is founcP^ to occur at about 110 K 
with a parasitic ferromagnetic moment whose remnant 
value is 0.025 /xs/Cr. The magnetic moment on each Cr 
atom was determined-^ to be about 2.6 hb and to point 
along the b two-fold axis but to rotate by 50° in the be 
plane as the temperature is lowered between 80 K and 
60 K. 

Kim et alW (KLVC) and Murakami et a/P (MFLL) 
have grown films of BCO with different results. KLVC 
grew films between 200 and 1400 nm thick on SrTiOa 
(STO) substrates with 8 nm thick epitaxial SrRu03 lay- 
ers as bottom electrodes. Their x-ray diffraction stud- 
ies revealed a tetragonal lattice with ao = cq = 3.888A 
and 60 = 3.902A. They find the Curie temperature T c = 
140 K and a weak ferromagnetic moment (FM) of 0.012 
/is/(f.u.) obtained from a sample cooled down to a very 
low temperature in a 2 kOe in-plane field. This weak 
FM can be compared with bulk crystal valued of 0.025 
/ie/(f.u.) obtained from a hysteresis curve. They also 
measured P(E), the electric polarization at 15 K as a 
function of a 1 kHz applied electric field in the direction 
normal to the film, and obtained the doublely looped 
hysteresis curve characteristic of antiferroelectrics. At E 
= ±700 kV/cm, where the P hysteresis loops closed, P 
had the rather small value of ±12/iC/cm 2 . The magnetic 
and electric polarization data reported by KLVC were on 
films 750 nm thick; which corresponds to a thickness of 
1922 single f.u. unit cells. MFLL grew BCO films on 
LaA10 3 (also on STO and NdGa0 3 but for which no 
data was reported). They obtained a triclinic structure 
with a = c = 3.9A, b = 3.88A, a = 7 = 90.6°, and 
ft = 89.1° which is only a small distortion away from 
being tetragonal and in fairly good agreement with the 
original room temperature ground state measurement. 
They obtained a weak FM of 0.05 ^s/(f.u.) from their 
strongly hysteric curve. Using piezoelectric force mi- 
croscopy MFLL find their film to be ferroelectric. This 
probably is consistent with their triclinic structure but 
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they could also be observing a surface effect. Except for 
a TEM image of a 200 nm BCO film, they did not report 
any film thicknesses. Recently, David et aZ.E2l reported 
that, using X-ray diffraction and transmission electron 
diffraction at 300 K, they found three phases in BCO 
films grown on SRO. The two they were able to identify 
were C2/c and Pnma. This is the first independent ex- 
perimental confirmation that Pnma BCO can exist below 
430 K. 

There exist four calculations for BCO of which we are 
aware, none of which contained the SOC necessary to 
obtain the parasitic ferromagnetism. Hill et aZp^ did the 
calculation for a hypothetical cubic perovskite structure. 
They found a soft phonon mode from which they con- 
cluded that BCO should be an A type antifcrrodistortive 
or antiferroelectric, i.e. ferroelectric planes with antifer- 
roelectric coupling between planes. They also predicted 
it was a G type antiferromagnet, i.e. every Cr atom 
surrounded by six Cr atoms with opposite spin. Xu et 
again found G type antiferromagnetism assuming 
the C2/c structure and calculated a Neel temperature 
more than two times larger than the experimental value. 
Baettig et al^ calculated the ferroelectric polarization 
and exchange coupling of BCO in the BiFe03 rhombohe- 
dral structure which has not yet been observed in BCO. 
Ray et aiP^ calculated stress in the prototypical cubic 
structure for different magnetic structures and studied 
the effects of magnetic ordering on the phonon frequen- 
cies in cubic perovskite structure. 

Our goal in this paper is to confirm the structure of a 
thick epitaxial film and to determine its properties. We 
do this by keeping the lattice constants fixed at those 
of the epitaxial film and calculating the positions of the 
atoms in the unit cell, finding the same Pnma structure 
as that when we allow the lattice constants to relax to 
their equilibrium values, and the same as the experimen- 
tal high temperature phase. The equilibrium and epi- 
taxial crystals are found to be antiferroelectric and an- 
tifcrromagnetic, but have spin-orbit induced spin tilting, 
the tilting direction depending upon a 2-fold axis along 
which the antiferromagnetic spins are originally aligned, 
but the direction of the D-M D vector is found to be 
independent of the original alignment. 



II. COMPUTATIONAL DETAILS 

Because the KLVC magnetic and electric measure- 
ments were conducted on films 1922 single f.u. unit 
cells thick, we feel confident in treating these films as 
bulk crystals. Calculations on the film were performed 
using the measured lattice constants, while determining 
the experimentally undetermined atomic positions within 
the unit cell. (By "film" we always mean relaxed atomic 
positions but unrelaxed experimental lattice constants.) 
The calculation was then repeated with both the lattice 
constants and atomic positions relaxed simultaneously to 
determine the Pnma ground state. For both calculations 



TABLE I. AE, the fully relaxed bulk crystal energy minus the 
film energy, in eV, calculated with four xc functionals. 





AE(No SOC) 


AE(SOC) 


GGA 


-0.04864 


-0.03966 


GGA+C/ 


-0.12025 


-0.11986 


LDA 


-0.52597 


-0.50464 


LDA+U 


-0.32089 


-0.30162 



simulated annealing with no symmetry constraints was 
applied. 

These calculations were performed using the projected 
augmented wave method^ as implemented in the VASP 
code^. For consistency with the International Tables^, 
we take a and c to lie in-plane, while b is normal to it. 
The four f.u. film lattice constants are a = c = \/2cq and 
b = 2&o with a and c rotated by 45°. The local density 
approximation (LDA) for exchange and correlation(xc) 
was treated with the Perdew-Zunger parameterization^ 
with and without a HubbarcP^ U- J of 2.2 eV. The gen- 
eralized gradient approximation (GGA) in the Perdew- 
Wang 91 forrrPS also with and without U-J — 2.2 eV was 
used as well. The outer core Bi 5d and Cr 3p electrons 
were treated as valence electrons. The energies with the 
SOC included were calculated self consistently using the 
atomic positions and lattice constants determined with- 
out the SOC. The atomic positions were iterated on until 
the largest force on any atom was 0.005 eV/A, resulting 
in a 10~ 6 eV convergence of the total energy. (All con- 
vergence tests were done for the GGA case.) All energies 
and lattice constants are per four f.u. (20 atom) unit cell. 
We used a 7 x 4 x 7 Brillouin zone sampling throughout 
except in one case where al0x6xl0 sampling revealed 
a convergence of 0.02 meV. A 500 eV cutoff in the plane 
wave expansion was used which resulted in the rather 
poor convergence of 70 meV in the total energy when 
compared with a 700 eV cutoff but a satisfactory con- 
vergence of 1 meV (2%) in AE, the bulk crystal energy 
minus the film energy. 

We do not discuss the bulk low temperature ground 
state because whether it is triclinic or monoclinic seems 
to depend on impurity content, imperfections, and do- 
main formation, which is beyond the scope of this paper. 



III. RESULTS AND DISCUSSIONS 

In table I we display AE, calculated using the GGA, 
GGA+U, LDA, and LDA+C/xc energy functionals. The 
GGA yielded the smallest energy difference, 40 meV (49 
meV without SOC), which is quite small for a 20 atom 
unit cell. Although small, this might be enough to desta- 
bilize a 1400 nm thick film, so, probably, actual value 
should be even smaller since the strained film has been 
observed. FIG.l. is the GGA film density of states, show- 
ing an energy gap of about 1 eV. All xc functionals re- 
sulted in semiconductors for both the film and relaxed 
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TABLE II. High temperature experimental and exper- 

imental filrrP^ lattice constants (in A) and volumes (in A 3 ), 
compared with four calculated results. 





a 


b 


c 


Volume 


Bulk 


5.5427 


7.7524 


5.4255 


233.130 


Film 


5.498 


7.804 


5.498 


235.939 


GGA 


5.612 


7.787 


5.442 


237.783 


GGA+U 


5.645 


7.816 


5.461 


240.931 


LDA 


5.455 


7.612 


5.344 


221.880 


LBA+U 


5.497 


7.647 


5.357 


225.194 



crystal. The lowest peak which contains 24 electrons per 
4 f.u. unit cell is almost pure oxygen s. The peak at -11 
eV, containing 8 electrons, mainly comes from Bi s. The 
isolated peak at the top of the valence band contains 12 
electrons. They are mainly Cr d electrons, which have 
mixed t2 9 and e g symmetry, and negligible O p ones. The 
72 electron band below it contains O p hybridized with 
Cr s, p, and d and with Bi s and p. The DOS includes 
SOC but the projections were obtained without it. 
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FIG. 1. Total GGA DOS(with spin orbit coupling) 



In table II we list the experimental bulk crystal, the 
film, and the calculated lattice constants, as well as 
the unit cell volumes. The experimental bulk lattice 
constants^ were determined at 470K where the ground 
state is orthorhombic with Prima symmetry and differ 
slightly from those at 500tP and 490KP. Note that a = 
c in the film, but not when the lattice constants are re- 
laxed, either theoretically or experimentally. The GGA 
relaxed volume however, is only 0.78% larger than the 
film's volume. The relaxed bulk crystal volume is larger 
than the bulk experimental volume by 2%, but if one 
tries to factor in what the experimental value would be 
if the bulk Pnma crystal structure were stable down to 
K, the error is more like 4% or 5%. This is typical 
for density functional calculations where errors of 2% in 
lattice constants are common. Note that with an appro- 
priate average of the GGA and LDA potentials, we could 



fit either the film or the experimental bulk volumes, get 
closer to the experimental lattice constants, and, thus, 
reduce AE and improve the stability of the film. 





a: 



r 



FIG. 2. (color online). Calculated structure of the BCO film 
showing the tilt of the O octahedra and the Bi displacements. 
In (a) one can see the displacement of the middle Bi plane 
with respect to the top and bottom Bi Planes. In (b) the same 
displacements are seen from a top view. Note this lattice has 
been translated by -b/4 relative to the coordinates used in 
table III. 



FIG. 2. displays the positions of the Bi atoms and the 
tilting of the oxygen octahedra around the Cr atoms in 
the film. Table III compares NPD measurements^ of the 
atomic positions with GGA calculations for the film and 
bulk crystal cases. Other xc functionals yield results sim- 
ilar to GGA. First note that, in spite of the fact that the 
film has a tetragonal lattice, both cases have the same 
Pnma space group as the high temperature experimen- 
tal phase. The Cr atoms are located at inversion centers. 
One can see that the Bi within a plane are displaced from 
their positions in an ideal perovskite lattice (in the GGA 
film, with similar results for other cases) by (0.04576a ± 
0.00874c) in one plane and by - (0.04576a ± 0.00874c) in 
the other, i.e. the interplanar coupling is antiferroelectric 
but the intraplanar a component is ferroelectric while the 
c component is antiferroelectric. Thus BiCrOa is a mixed 
A and G type antiferroelectric. The Bi displacement to- 
gether with the rotation of the oxygen octahedra are the 
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TABLE III. GGA calculated film and fully relaxed struc- 
ture parameters together with experimental high temperature 
data^. See Ref.19 for Wyckoff notation and other atomic po- 
sitions. 





atom 


Wyckoff 


a 


b 


c 


Experiment 


Ol 


8d 


0.2948 


0.0390 


0.7015 




02 


4c 


0.4793 


1/4 


0.0811 




Bi 


4c 


0.0432 


1/4 


-0.0051 




Cr 


4b 


1/2 


1/2 





Film(GGA) 


01 


8d 


0.29356 


0.04390 


0.70571 




02 


4c 


0.48092 


1/4 


0.08687 




Bi 


4c 


0.04576 


1/4 


-0.00874 




Cr 


4b 


1/2 


1/2 





Fully relaxed 


01 


8d 


0.29546 


0.04291 


0.70153 


(GGA) 


02 


4c 


0.47696 


1/4 


0.08755 




Bi 


4c 


0.05141 


1/4 


-0.00927 




Cr 


4b 


1/2 


1/2 






TABLE IV. Energy (in meV) of BCO film for different spin 
orientations relative to the 6 axis, i.e. to the (010) direc- 
tion, calculated in the GGA and GGA+U. 

oriention GGA GGA+U 

Oil) 

100 0.162 0.039 
001 0.160 0.003 

101 0.153 0.006 
011 0.077 0.005 
110 0.072 0.004 



origin of the antiferroelectroic nature of Prima BCO. 

We found that for all four xc functionals, the magnetic 
ground state is, as expected, G type. We found the C 
type, A type, and ferromagnetic states are 35.93, 153.76, 
208.84 meV above the G type, respectively. These results 
are GGA without the spin-orbit interaction. The SOC 
was introduced with all the Cr spins initially aligned in 
one of several directions, some of whose energy relative 
to the b alignment are displayed in table IV. Although 
the SOC energy is 2.02 eV, the anisotropy energy dif- 
ferences are small and the spins remained in their origi- 
nal orientations, except for the spin-orbit induced devi- 
ations displayed in table V. When the SOC was turned 
off but noncollinearity still allowed, the spins returned 
to their original collinear antiferromagnetic orientations. 
The magnetic moment of 2.64 \ib on each Cr agrees well 
with the value of 2.6 /j,b deducecP from NPD data, albeit 
for a different crystal structure. We note that for the b 
orientation, a and c components are induced on each Cr 
atom but the a components are antiferromagnetic and 
therefore cancel, leaving a SOC induced ferromagnetism 
of 0.073 fiB per Cr in the -c direction. Similarly, the c 
orientation has 0.079 /is induced ferromagnetism along 
-b with the a components again cancelling. But for the a 
orientation the b and c components both cancel and the 
crystal remains antiferromagnetic. The fact that in no 
case is there any ferromagnetism along a complicates the 



TABLE V. Local spin and orbital moment components for dif- 
ferent initial antiferromagnetic spin orientations in the GGA. 

orientation spin orbit moment 

a b c a b c 



Crl 


0.013 


-2.637 


-0.073 


-0.001 


0.042 


0.007 


Cr2 


-0.013 


2.637 


-0.073 


0.001 


-0.042 


0.007 


Gr3 


0.013 


2.637 


-0.073 


-0.001 


-0.042 


0.007 


Cr4 


-0.013 


-2.637 


-0.073 


0.001 


0.042 


0.007 


Total 








-0.292 








0.028 


Crl 


2.637 


0.028 


-0.058 


-0.044 


-0.004 


-0.005 


Cr2 


-2.637 


-0.028 


-0.058 


0.044 


0.004 


-0.005 


Cr3 


-2.637 


0.028 


0.058 


0.044 


-0.004 


-0.005 


Cr4 


2.637 


-0.028 


0.058 


-0.044 


0.004 


-0.005 


Total 




















Crl 


0.043 


-0.079 


2.637 





0.009 


-0.040 


Cr2 


0.043 


-0.079 


-2.637 





0.009 


0.040 


Cr3 


-0.043 


-0.079 


-2.637 





0.009 


0.040 


Cr4 


-0.043 


-0.079 


2.637 





0.009 


-0.040 


Total 





-0.316 








0.036 






interpretation of the 0.012 [Ib measured valuers which 
was obtained by applying a magnetic field along one of 
the in-plane <100> directions of the substrate. It is 
likely, that to reduce the strain, there are domains in 
which a is in the <100> direction and others in which it 
is in the <001> direction which would make the exper- 
imental film value consistent with the theoretical value 
averaging over different domains though the theoretical 
value is a little bigger. Additionally, we expect reduced 
energy denominators, caused by GGA density functional 
gap narrowing, to enhance the SOC and therefore en- 
hance the calculated SOC induced ferromagnetism, thus 
making the disagreement between the calculated and ex- 
perimental ferromagnetism understandable. In fact the 
GGA+ U energy gap is about twice^ that of the GGA, 
resulting in initially b(c) antiferromagnetic spins having 
0.049 (0.053) [is induced ferromagnetism in the c(b) di- 
rection. 

We are now in a position to see how the D-M theory 
applies to the current case. We note that the a coordi- 
nate is unique in that no weak ferromagnetism is induced 
along it if the antiferromagnetic spins oriented along ei- 
ther of the other axes, or when a is the antiferromagnetic 
axis, there is no weak ferromagnetism along the other 
two axes. This uniqueness is a consequence of the dif- 
ferent non primitive translations^ associated with the 
two- fold rotations about the different axes. The transla- 
tions belonging to the b and c rotations interchange (G 
type) spin sublattices whereas the translation belonging 
to the a rotation does not. The unitary Prima symmetry 
is preserved only when the antiferromagnetic polariza- 
tion is oriented along a. Because each component M a of 
the ferromagnetic vector M is rotated into -M a by one 
of the perpendicular two-fold rotations (the nonprimi- 
tive translations have no effect on M), ferromagnetism 
is forbidden (unless there is a broken symmetry). When 
the antiferromagnetic polarization is oriented along the b 
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axis, the unitary group is reduced to E, I, <r c [a/2fi,c/2], 
and U c [a/2,0, c/2]. These four operations plus r, the 
time reversal operator, multiplying the remaining four 
reflection and rotation operators, cr&[0, 6/2,0], Uf, [0,6/2,0], 
a a [a/2,b/2,c/2], and U a [a/2,b/2,c/2] form the antiferro- 
magnetic antiunitary group of the crystal. All eight oper- 
ations leave the c (but not o) component of the spins' sign 
unchanged. If the antiferromagnetic polarization is along 
c, interchange the 6 and c subscripts to obtain the mag- 
netic group. Then all eight operations leave the b (but 
not a) spin component unchanged. Therefore, all oper- 
ations of the magnetic group are consistent with weak 
ferromagnetism along the b or c axis, respectively, if the 
antiferromagnetic spin orientation is along the c or b axis. 
One can also show that the antiferromagnetic spin tilts 
belong to the magnetic groupP^i. If a D-M DSi x S2 term 
is present in the Hamiltonian, it is easy to see that D 
must be polarized in the a direction, irrespective of the 
orientation of the spins. 

In summary, we have performed density functional cal- 
culations for a BCO film with lattice constants deter- 
mined by the STO substrate on which it was epitaxially 
grown as well as for a bulk crystal with relaxed lattice 
constants, and found in both cases that they have the 
Pnma space group of the high temperature phase, indi- 
cating that the high temperature phase is stable (when 
grown epitaxially on STO) down to the 10 K to which 
the film was taken. We obtained a value for the mixed A 
and G type antiferroelectric Bi displacements; We found 
G type antiferromagnetism and if the antiferromagnetic 
spins were aligned along the 6(c) axis a weak ferromag- 
netic moment was spin-orbit induced along the c(6) axis 
and if they were aligned along the a axis there was no in- 
duced ferromagnetism. The uniqueness of the a axis was 
attributed to the nonprimitivc translation accompanying 
its 2-fold rotation which does not interchange the spin 
sublattices. It was found that SOC induced antiferro- 
magnetic spin tilting, consistent with the magnetic space 
group, always occurred where ferromagnetic tilting did 
not. We showed how this could have been predicted from 
a theoretical analysis of the Pnma group, and that D in 
the D-M Hamiltonian must point along the a axis irre- 
spective of the antiferromagnetic spin polarization, Note, 
however, that does not account for the antiferromagnetic 
spin tilting. 
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